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Charge Neutralization in an Ionic Rocket 


Cue iClouten andl owe deka ens 


bstract 


An ionic rocket model is considered in which ions are 
ected from an ideal plane diode. Charge neutralization of 
e ion beam is accomplished by means of electrons injected 
ound the ion beam. The potential within the ion beam in the 
esence of electrons is studied as an electrostatic problem. 
<ial and radial solutions are given for various charge distri- 
tions. It is found that the potential at the centre of the beam 
pends primarily on the charge per unit length unless com- 
ete neutralization is accomplished very rapidly. 

The power and thrust of an ionic rocket is studied in terms 
the parameters of an ideal plane diode. Considerations of 
e potential depression within the diode indicate that the 
ywer and thrust deteriorate very rapidly for diode spacings 
ss than one-half the radius of the beam. For these same diode 
acings, the reduction in the ion diode efficiency due to the 
sctrostatic force of the partially neutralized beam becomes 
ry important. 

Considerations are given for the optimum choice of diode 
wameters. A numerical example illustrates the relative im- 
tance of the various factors. The process of charge mixing 
discussed in view of the approximations made. 


1troduction 


It is well known that the useful part of the exhaust 
an ionic (i.e. ion diode) rocket must contain equal 
nounts of positive and negative charge and that any 
cess will be electrostatically attracted back to the 
cket and thus contribute no thrust. If a planar diode 
infinite cross section is considered as the diode 
celerator, it seems that intimate charge mixing should 
ke place within a distance which is small compared 
th the electrode spacing. For practical current 
nsities and voltages this spacing will be very small 
d mixing within this distance very difficult if not 
possible. 
This limitation is however only severe in an infinite 
de. In a diode of finite area the allowable distance 
charge separation is increased as the ratio of diode 
a to square of the spacing is decreased. The actual 
metry of the ion and neutralizing electron beams 
1 also be more readily treated in a model with a 
ite area. 
t will be seen that the space charge effects in the 
de are of the same order as the exhaust charge- 
aration effects for an electron beam with area of the 
e order as the ion beam area; both are reduced by 
ucing the area for a given spacing. The ideal thrust, 
a charge effects, for a given diode spacing and 
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voltage varies directly as the diode area, but reductions 
from the ideal thrust due to space charge effects are 
also more noticeable as the area is increased. In other 
words, a very large area diode with severe charge limita- 
tions can, with advantage, be subdivided into smaller 
diodes until the point of diminishing returns is reached, 
where further subdivision produces negligible improve- 
ment in thrust for the added complexity of the struc- 
ture. 

The charge-separation effect in the exhaust will be 
examined by considering an ionic rocket in which ions 
are ejected from an ideal plane diode of circular cross- 
section. It will be assumed that the charge neutraliza- 
tion is accomplished by means of electrons injected 
around the ion beam. Because of the electric field 
across the ion beam the electrons will oscillate radially 
through the beam with the result that they will form 
a cylinder of negative charge concentric with the ion 
beam. The potential within the ion beam in the presence 
of electrons will be treated as the electrostatic problem 
of two concentric cylinders of opposite charge. 

The power and thrust of an ionic rocket will then be 
studied in terms of a parameter (electrode spacing/ 
beam radius) of an ideal plane diode. Now both the 
potential depression (charge effects) within the ionic 
diode and the ejection-plane image force due to in- 
complete charge mixing outside the diode in the exhaust 
will be considered. 

This analysis leads to the conclusion that the elec- 
trode spacing should be in the range of 0.5 to 1 times 
the beam radius. It is also deduced that complete 
charge mixing and/or recombination further than a 
distance of about ten beam radii has little effect on the 
system. From a typical example it is found that com- 
plete neutralization will be negligible within such a 
distance. 

Since electrical and mechanical quantities are in- 
volved in several expressions, it was found convenient 
to use the M.K.S. system throughout the paper. 


I. Electrostatic Solutions 


1. Axial solution for two concentric cylinders of opposite 
charge in front of an infinite conducting plane 


The model to be studied is represented by two 
concentric cylinders of opposite charge placed in 
front of an infinite conducting plane. To obtain the 
solution for the potential along the axis of the cylinders 
the conducting plane is replaced by two concentric 
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image cylinders as shown in Fig. 1. Since the potential 
is a scalar quantity, the final solution may be obtained 
by proper summation of the potential due to each 
cylinder. 


a) Uniform charge distribution 


Single cylinder. The potential at a point 2 on the 
axis of a finite cylinder of positive charge is: 


Vv. = 1 ce Qv 2ar dr dx Ge) 


hea (x — a)® + 27} 


where gq, is the volume charge density 
ro is the radius of the cylinder 
Lis the length of the cylinder 
€) is the permittivity of free space. 
After integration, if one defines 
ga = Qo Tre, the charge per unit length 
A= 1/ fy ANOEX “= Boy 7a 
the expression for the potential along the axis becomes:T 


7 A detailed derivation of this solution has been shown in 
the Research Report No. 7-811-1, Aug. 1959, RCA Victor Co. 
Ltd., Montreal, Canada. 


Vee a ANS Xr a ead 


chp Foe com pe ah ce 
seisinhon(Al-—-#.X)¢-tsinhen( X)| 


Ion and electron cylinders and their images. The 
axial solution for the potential of two concentric 
cylinders of opposite charge in front of an infinite 
conducting plane is obtained by superposition of four 
equations of the type of Eqn. (2). 


Ve = Vas F View + Vio. ++ Va, (3) 


where V, is the total solution (See footnotef) 

V., is the solution for the cylinder of positive 
charge (72) 

V., 18 the solution for the cylinder of negative 
charge (r;) 

V.,18 the solution for the image cylinder of 
negative charge (r2) 

V., is the solution for the image cylinder of 
positive charge (7) 


=— qa ies a7 Ne + wry? 14 
Vo, = gt | (a = wONCA = XO AP XE? 4 


ae (Arex Th ee eee (3a) 


re 


+ sinh’ (A — X) + sinh7 2x) | 


A Trey 
+ RoX[RoX" + 1} — Ro'(A — X)|A— X| 
— Ry X |X| + sinh Roi(A — X) (3b) 


V, =o Ex — X)[Ro(A — X) + 1} 


+ sinh” (FX) | 
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Ve, = 
: Atreg 

+ (A+ X(A+ XP +17 + X|X| 

— (A+.X)|A+ X|- sinh™ (X) 


=a | - xix? + 1 


( 
+ sinh (A + x)| 


Ve | - Ry X{Re'X? + 1} 


+ R(A + X)[Ro(A + X)* + IP +. ReX | X | cs 
— Re(A + X) (A + X= sinh Gros 
| 


+ ginh, fata x)| 


In these expressions, the length (A) of the cylinder 
and the distance (X) along the axis have been not 
malized to the radius (r2) of the cylinder of positiy 
charge. The ratio of the radii of the two concentri 
cylinders is expressed as Ro = 11/12. a 

The solution has been evaluated for A = l/re : 
0.3, 3 and 30 and for Ry = 1.5 and is shown on a log-lo: 
scale in Fig. 1. It is noted that for A = 30 (long cylinder 
the conducting plane has little effect on the potentie 
throughout the cylinders. The potential attains th 
same value as if the cylinders were infinitely long 
However the situation is quite different if the cylinder 
are short (A = 0.3) that is when complete neutraliza 
tion of the ion beam is rapidly accomplished. 


b) Linearly decreasing charge density along the azi 
To allow for gradual mixing of the charges th 
problem of concentric cylinders of charges has bee 
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DISTANCE ALONG THE AXIS ( X) 


Fria. 1. Potential variation along the axis of two concentr 


cylinders of opposite charge in front of an infinite conductit 
plane. 


udied for the case where the charge density is uniform 
ver the cross section but linearly decreasing along the 
sis. The solution is obtained as above (Eqn. 3) by the 
ethod of images. (See Footnotet.) 

A plot of the potential variation along the axis for 
le case A = 30 and Ry = 7/r2 = 1.5 is shown in Fig. 
in comparison with the equivalent solution for the 
ise of uniform charge density throughout the cylinders. 
Ithough the potential drops more rapidly when the 
1arge density decreases linearly it is seen that, for 
= 30, the maximum value of the potential reaches 
/% of that for an infinitely long cylinder. Therefore, 
the neutralization occurs in a distance greater than 
) times the radius of the beam, the theory for uniform 
1arge distribution along the axis is adequate. As will 
> shown later, practical recombination or mixing 
stances will give values for A much greater than 30, 
y curves for A = 0.3 and 3 were not calculated. 

_ Radial solution for two infinitely long concentric 

cylinders of charge 


Let us consider two concentric infinite cylinders of 
oposite charge. The small cylinder contains positive 
1arges uniformly distributed throughout the volume. 
1 the large cylinder, the negative charge density is 
ssumed to increase radially as the square of the radius. 
he solution for the potential across the two cylinders 
expressed in terms of two parameters: the relative 
dii of the cylinders (Ro) and the ratio (vy) of the 
avative charge densities at the centre and the edge 
' the large cylinder. Let the charge densities be 
ee Fig. 3): 
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ae 
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LINEARLY DECREASING 
CHARGE DENSITY 


ie} 4 8 v2 16 20 


UNIFORM 
CHARGE DENSITY 


and 


de, = 4 


Up} 
where q,, is the negative charge density at r = O and 
a is a parameter defined in terms of the ratio (y) of 


the negative charge density at the centre and the edge 
of the cylinder. 


YS qvo/ Tor, (5) 


Equation (4) becomes: 
r 2 
= NG) | 
ry 


Go, = —~ Qo E +( 
7 


dv) 18 defined by assuming the same charge per 
unit length for both cylinders: 


a 


TL 
= —2n | Qo, (r)r dr 
0 


which gives upon integration: 


qa 1 

ee my [1 + 1/2(1/y — 1)] - 
The electric field across the two concentric cylinders 
is obtained by applying Gauss’s theorem. The potential 
is then derived by integration of the electric field along 
the radius of the cylinders. The constants arising from 
the integration are determined with the help of the 

following boundary conditions: 


Vax) =O 

ote Mi Aer) 

Defining R = r/rz and Ry = 1/re (i.e. normalized as 
before to the ion beam radius) one obtains the follow- 
ing expressions for the potential across the cylinders. 
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me ot 
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1G. 2. Potential variation on the axis of two concentric cylinders of opposite charge in front of an infinite conducting plane. 


iform and linearly decreasing charge densities. 
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For Orne ii 
V(R) = ae Inthe +p, (l=? he?) 

(= RRs) 

+ ig0/y — DI ny 

_ O/y- Da - pay 
AE aL) 
For erecta 
_  (— R/Re) 

VR) = Be] 2m Rel R) — Gage DT (7) 


_ G/fy = DU = Ee) 
afl + 61/7 — DI 


The variation of potential across the two cylinders 
has been plotted for two limiting cases y = O and 
y = 1 (Fig. 3). The case y = O represents the condi- 
tions for which most of the electrons are concentrated 
near the outer radius, while for y 1 the electrons 
are assumed to be distributed uniformly across the 
cylinder. 

The mean value of the potential over the cross sec- 
tion of the positive cylinder may be calculated by 
means of the following relation: 


ve a V(R)2aR dR (8) 


For Ry = 1.5, with y = 0, V is 0.644 Vo and with 
y = 1, V is 0.657 Vo where Vo is the potential on the 
axis. 
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II. Power and Thrust of the Ionic Rocket 


Consideration will be given to the effective powe 
and thrust of an ionic rocket first in terms of the pi 
tential depression due to the space charge in the diod 
and then in terms of the image force due to incomplet 
charge neutralization of the ion beam. 


1. Potential depression in an ideal plane diode 


Due to the space charge in the ionic diode a radig 
as well as an axial potential depression exists in th 
beam. The axial depression limits the actual curren 
drawn by the diode (space charge limited flow). Th 
radial potential depression reduces the energy of th 
ions in the centre of the beam. H 

In order to estimate the effective power and thrus 
of the diode in terms of the potential depression, 
will be made of the theory given by Smith and nal 
man (1) and Pierce (2) for the case of an homogeneou 
beam of circular cross section. In the derivation of th 
theory the current density is assumed to be uniforn 
across the beam while the charge density and the volta 
vary. In fact, the current at the centre of the beam i 
smaller than that at the edge. However, since th 
reduction in the current at the centre of the beam wi 
induce a reduction of the potential depression acros 
the beam, it is believed that the above assumption i 
not unreasonable for the calculations of the effective 
power and thrust. Reduced emission at the centre 0 
the emitting plane and reduced voltage in the aperture 
(lens effect) will not be considered here. 

Using Pierce’s notation, a potential ratio is definec 


4 


Z| 


Fig. 3. Potential across two concentric infinite c 
density (y = 0). Uniform negative charge density ( 


ylinders 
Y S10 
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ie} 


NORMALIZED RADIAL DISTANCE , R 


of opposite charge. Maximum depression of the negative charg 


the ratio of the voltage at the beam edge to the 
ltage in the centre of the beam: ~ 


¢=V/V, 


rere Vj, is the voltage at the centre of the beam 

V is the edge voltage 
erce (2) has plotted ¢ as a function of a reduced 
dius p defined in the following way: 


‘ 1 raNe | WN 
— meo|2q/M |? (<) (;) 


here r is the beam radius. 
a is the radius of the conducting cylinder 
containing the beam 
In the present case ‘‘a”’ is the radius of 
the exist aperture, so V = Vo the electrode 
voltage. 
€« is the permittivity of free space 
A I ~ is the total current 
q/M is the charge to mass ratio of the ions; for 
singly charged ions, q is the electronic 
charge. 
ot M = 1836 mA, where m is the mass of the electron 
id A is the atomic weight/number electron charge 
id r = a, i.e. the beam extends over the whole aperture. 
For the situation of interest in the present problem 
here the potential ratio ¢ is close to unity, it is 
sumed that V; V and the ion perveance can be 
ritten as: 


(9) 


fort 


Pi= I/Vo3!? ~~ 1/V3! (10) 


; may be expressed in terms of the electron perveance 
of the same diode with electron flow a function of the 
ructure only, 


Sooke (11) 


so that the dependence of P; on the ion mass is shown 
explicitly. For a plane diode of circular cross section, 
the electron perveance is: 


_ 7.35 X LO 


fg? DP 


(12) 


where D = d/r 
d is the diode spacing 
r is the radius of the exit aperture (ion beam 
radius). 
Again d is normalized to the ion beam radius. Using 
the relations (10), (11) and (12), equation (9) becomes 


2_ 735 X10° 1 _ 445 
‘ (2e/m)? De 


D 
It is more convenient to define the potential depres- 
sion as the ratio of the potential difference between 
the edge and the centre of the beam to the potential 
at the edge: 


(13) 


Vi —il 

: lies , e 

The potential depression (1-¢~!) has been plotted 
(Fig. 4) as function of the normalized diode spacing 
(D) making use of Eqn. (13) and the curve of ¢ vs. p 
given by Pierce (2): 

Effective Power Ratio. The effective power of the 
diode can be defined as: 


Wa = 1V.= 


Potential depression = 


PVee 
428A} 


where J is the total current 
V.is the potential of the beam averaged over 
the area of the exit aperture. (Strictly 
speaking V, should be the square of the 
average of V1 over the aperture, but since 
V, is close to Vo for all useful cases, the 
error in this average V is small.) 


i eal eee 
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Fra. 4. Potential depression in the exit aperture of the diode as function of the normalized diode spacing 
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The effective power ratio is defined as the ratio of the 
power taking into account the potential depression 
across the beam to the power if the potential were 
uniform across the beam. 


Effective power ratio 


Wee 
due to diode effects = nwa = 1 2 V, 


Let the average of the potential depression over the 
cross-section be a fraction (a) of the maximum depres- 
sion. Then: 


Va =e Vo = a(1 = o)Vo (16) 
and 
wa = Feigiy ee (17) 


The percentage effective power ratio has been plotted 
in Fig. 5 as function of the diode parameter D using 
Eqn. (7) and the curve (1 — ¢~) vs. D of Fig. 4. For 
this curve a was chosen as 0.657 (see Eqn. 8). 

Effective Thrust Ratio. The effective thrust of the 
diode can be written as: 


eMC, 
g 
where the thrust is expressed in. Kilograms 
Ce is the ion exhaust velocity averaged over the beam 
M is the mass flow rate 
g is the acceleration due to gravity 
In terms of current and voltage, 
Lf 2M gee 
Rg = >| —)}) 1IV7 == — 18 
é (=) g 478A? Mie) 


Is the effective thrust ratio (nr,) due to diode effects 
if defined as the ratio of the thrust taking into account 


Fy 


EFFECTIVE RATIO, % 


the potential depression in the beam to the thru 
assuming uniform potential across the beam, the 


Fa Ve - 
py ex Ee 1 
Fa Fo G) ( 
Defining V, as in Eqn. (16), 7” becomes: 

NF a — {1 = a(1 == aoe ( 


A curve giving the percentage effective thrust ratio ; 
function of the normalized diode spacing is shown 
Fig... we ; 

Actual Thrust and Power. The ideal thrust of — 
diode with uniform potential across the beam can b 
expressed as: ; 


Sant ; 
Pree ; ( x 42.847IV 3 (21 


and in terms of the electron perveance for a plane diod 
(using Eqns. 11 and 12) 


Fy = 735 X 10° C2) Vor 
g\e D» 
oy (22 
= 252 X10 — = Kgm 


Curves representing the thrust (/o) as function of J 
are shown in Fig. 6 for accelerating voltages of 5, 1¢ 
and 20 KV. The effective thrust (Fz) of the diode ca 
also be calculated by means of equations (19), (20) an 
(22). Fa is also shown in Fig. 6 in comparison with th 
curves for Fy vs. D. The power of the diode assumin 
uniform potential across the beam can be expressed as 
(using Eqns. 11 and 12) 


V 5/2 
r 0 
A?D? 


Combining equations (15), (17) and (23) the effectiv 


W. = IV = PV’? = 172 X 10° Watts (23 
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Fra. 5. Effective power at ‘ust ratios of ionie di ; i i i 
> I or and thrust ratios of the ionie diode as function of the normalized diode spacing 
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Fig. 7. Sketch of the potential variation in an ionic (diode) 
rocket. 


power may be obtained as function of D; curves similar 
to the ones for the thrust of the ion diode are then 
obtained. 

2. Effect of the image force due to incomplete charge mixing 


To estimate the efficiency of the diode as a thrust 
device in free-space, it is necessary to take into account 


4 


Fp (IDEAL THRUST) 
—— —F,, (EFFECTIVE THRUST) 


6 1:8 


NORMALIZED DIODE SPACING, D 
Fia. 6. Ideal and effective thrust of the ionic diode as function of the normalized diode spacing (for cesium ions) 


the reduction of thrust due to the electrostatic force 
exerted on the diode by the partially neutralized ion 
beam. In view of the results obtained above, the overall 
distribution of potential in the system can be qual- 
itatively represented by sketches of Fig. 7. In addition 
to the loss of energy within the diode due to the space 
charge, the ions suffer another loss of energy outside 
the diode due to the electrostatic force between the 
incompletely neutralized ion beam and the diode itself. 
Since the problem is treated as an electrostatic one, the 
value of the potential in the ion beam outside the diode 
represents the energy necessary to keep the ions in the 
configuration assumed at the origin of the problem. 
The maximum height of the electrostatic ‘hill’? can 
be expressed in terms of the plane diode parameters. 
As a result of Eqn. (3) and Fig. 1, the maximum 
value of the potential function in the case of concentric 
cylinders of opposite charge can be expressed as: 


(24) 


is the ion charge per unit length before 
neutralization 
f(A) is the maximum value of the potential 
curves (Fig. 1) for a given neutralizing 
distance. 
Values of f(A) have been calculated for A = 30, 6, 3, 
1.5, 0.9 and 0.3. The resulting curve is shown in Fig. 8. 


where qa 


Since 


and 
- ry 3/2 
I = PV 0 ) 


for the zdeal plane diode, the expression for qa in terms 
of the plane diode normalized spacing becomes, using 
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Eqns. (11) and (12): 


—6 fm ; Vo 
= 7. ed pee 25 
qa = 735 X 10 @ ri (25) 
Substituting this value of gq, into Eqn. (24) gives 


Van 725 10 (ee f(A) 
pete ee eae es) Sikes ae 
Ate @ Dz pant Do ( 6) 


The effective exit power and thrust ratios (relative to 
the ideal diode) due to the image force of the partially 
neutralized beam may be expressed in a manner similar 
to that of the preceding section (Eqns. 15 and 19). 
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By means of the function f(A) (Fig. 8) and the valu 
of Vy (Eqn. 26), the effective thrust ratio nr, ha 
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Fra. 8. Variation of the maximum potential on the axis of the ion beam as function of the neutralizing distance (uniform charg 
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ffeen plotted in Fig. 9 has function of D for various 
values of A. 

_ As seen from Fig. 9 no improvement in the effective 
thrust ratio is expected unless the charges are mixed 
very quickly, i.e. for A smaller than 10. 

_ For velocities at which the ions are ejected from the 
diode, it seems very unlikely that the charges can be 
Mixed or recombined within such a short distance. 
For example a beam of 0.5 em. radius has to be mixed 
within 5 cm. (A = 10) from the exit aperture. For a 
6.8 KV beam of cesium ions this distance correspond 
to an unrealistically short mixing time of 0.15 usec. 

For practical purposes the model of the two long 
concentric cylinders of charge therefore represents an 
adequate approximation. 

It is seen from the curves of Fig. 9 that the un- 
desirable electrostatic action of the incompletely 
neutralized ion beam becomes negligible when the 
diode spacing is allowed to become large (low perveance, 
low space charge, low mass flow rate). However, from 
the results presented in Fig. 6 it is observed that the 
diode spacing has to be made as small as possible for 
obtaining useful mass flow and thrust from the diode. 
It appears that a diode perveance corresponding to a 
spacing D between 0.5 and 1.0 would give the best 
compromise. For values of D smaller than 0.5, the 
effective thrust ratio (Fig. 9) drops very rapidly, due 
to excessive charge effects. If the value of D is allowed 
to become larger than 1.0, the actual power and thrust 
would become too low for practical applications. 


Il. Practical Considerations 


A numerical example may illustrate the use of the 
formulae. The given parameters from Ref. (3) are: 


Ton Current 


Density: J) 10'ma/enm? = 10? Anip/m?2 
Ion Beam Radius: 7 = 0.75 cm 2 ei) mM 
Electron Injector 

Radius: fp == ).o°Cm = 0- x 10 m 

or hy = 1.5 
Applied Voltage: Vo = 6.8 KV 


Singly Charged 


Cesium Ion: A 133 


Potential at the centre of the ion beam. The potential 
at the centre of the ion beam (Rk = 0) with respect to 
the potential on the accelerating electrode may be 
calculated by inserting the above data into Eqn. (7). 

For the case y = 0, i.e. maximum depression of the 
electronic charge density, the potential drop is 925 
Volts. 

For the case y = 1, ie. uniform electronic charge 
density the potential depression is 572 Volts. 

Stuhlinger and Seitz (3) have studied the motion of 
an electron injected at a given distance from the axis of 
an infinitely long ion beam. They have estimated the 
energy of the electron on the axis of the beam by 
correcting for the average time spent by the electron 


within the ion beam. Using Stuhlinger and Seitz’s 
theory, the effective value of the potential drop at the 
centre of the ion beam for this case is (after correction 
of an error in Stuhlinger and Seitz’s formula) 820 Volts. 

It is found that by solving Eqn. (7) for R = 0, their 
value corresponds to a parabolic electron density with 
the electron density at the centre 0.185 times the edge 
density. 

Equivalent value of D. By means of equations (10), 
(11) and (12) an equivalent value of D may be calcu- 
lated for the numerical example given above. D is 
found to be 1.03 which represents the upper limit for 
the optimum choice of this parameter. 

Characteristics for a typical ionic rocket. For a diode 
spacing D of 1.03, the effective thrust of the diode is 
2 X 10-* Kem (See Fig. 6). The reduction of thrust 
and power due to the electrostatic action of the partially 
neutralized ion beam is only 3% and 6% respectively. 
For a rocket system weighing 1000 Kgm and having a 
specific weight of 10 Kgm/Kw (just practical at present 
without payload) it would be necessary to use 1000 
diodes of this type to obtain a thrust of 2 * 107 Kgm 
giving to the vehicle an initial acceleration of 2 X 10~*g. 

Imtial condition of charge mixing. In the above 
analysis it was assumed that the electron distribution 
through the beam was established immediately at the 
origin (diode accelerating electrode). The validity of 
this approximation may be tested by computing the 
time taken by an electron to reach the centre of the 
beam. This can be done using Stuhlinger and Seitz’s 
theory and a value of 6.95 * 10~° see is obtained for 
the set of data given above. 

For a velocity of 100 Km/sec (6.8 KV) along the 
axis of the beam, the cesium ions travel only 7 X 107? 
cm before the electrons reach the centre of the beam. 
Therefore the electron distribution across the ion beam 
takes place effectively at the exit aperture of the diode. 

Beam Spreading. In the above analysis, it has been 
assumed that the ions move parallel to the axis of the 
system. In fact, it is well known that due to the space 
charge forces the ion beam will spread out. The spread- 
ing of the ion beam will have the effect of reducing the 
thrust of the diode if the ions and electrons are not 
completely mixed prior to a distance for which the 
divergence of the beam is too large. 

Since the thrust will vary as the cosine of the angle 
of divergence, a large spreading of the ion beam may 
be tolerated before an appreciable reduction of thrust 
is obtained. For a beam divergence of 30° the thrust at 
the edge of the beam will be reduced by 13% of its 
value if no spreading were occurring. The overall 
reduction of thrust should be even smaller since the 
angle of divergence for the ions within the beam is less 
than for the ions at the edge; for a uniform distribution 
the overall decrease is 11%. Furthermore the beam 
spreading is likely to facilitate the mixing of electrons 
and ions, and consequently should tend to reduce the 
image force of the ion beam. 


o1 
a1 


The Journal of the Astronautical Sciences 


Recombination and mixing distance. It was found, in 
the above example, that the electron energy at the 
centre of the ion beam is in the range of 600 to 900 eV. 
Therefore at the beginning of the mixing process the 
beam will effectively behave like a hot plasma in which 
the electrons will lose their energy mostly through 
inelastic collisions with the ions. Dandurand and Holt, 
(4) have reported a value (ac, = 10-7 em’/ion-sec) for 
the recombination coefficient between electrons and ions 
in the afterglow of a cesium discharge (about 1 electron 
volt energy). The electrons in the plasma from an ionic 
rocket would be much faster (hundreds of electron 
volts) than those for which Dandurand and Holt ob- 
tained their value for acs. For the above numerical 
example, the ion density is 6.25 X 10° ions/em* and 
the recombination time is 1.6 milliseconds for slow 
electrons and much longer for fast electrons. A beam 
of 6.8 KV cesium ions will travel 0.16 Km, during that 
period, and recombination would have only just begun. 

It is not essential, however, that the ions and 
electrons be completely recombined since as soon as the 
ions and electrons are being mixed, the electrostatic 
potential within the beam will decrease. This situation 
is illustrated in Fig. 3 where the potential on the axis of 
the beam is plotted for various values of the ratio (Ro) 
of the radii of the electron and ion beams. This ratio 
Ry can be considered as a measure of the amount of 


mixing between ions and electrons. The curves of Fig. 
also indicate that there is a definite advantage o 
injecting the electrons as closely as possible to the io1 
beam. 


Conclusion 

As long as the overall neutrality of the ionic rocke 
exhaust is preserved, complete charge mixing and/o 
recombination further than a distance of ten bean 
radii has little effect on the system. Thus only charg 
separation near the rocket affects the performance. 

It has been found that an electrode spacing betwee! 
0.5 and 1 times the ion beam radius represents th 
optimum range of operation of an ideal plane diode a 
a thrust device. Above this range the mass flow rat 
and the thrust are too low for practical applications 
For a spacing below this range the space charge effect 
limit the performance of an ionic diode quite sharply. 

This work was supported by Astro-Electronic Prod 
uct Division of Radio Corporation of America 
Princeton, New Jersey. 
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phemeris Time 


T. P. Gabbard* 


Abstract 

EKphemeris Time is defined and its meaning briefly dis- 
cussed ; following this is a historical development of the events 
leading to the recognition of Ephemeris Time as a successor to 
Universal Time for use in precise work. The determination of 
Ephemeris Time utilizing lunar datum and atomic and/or 
molecular clocks is delineated, and some uses of the more 
uniform time measure are indicated. 


1. Introduction 


Man has long been concerned with time, for it has 
regulated his life from his first dawn. Time is the 
subject of this paper; and in particular, man’s latest 
achievement in denoting the uniform passage of this 
fleeting dimension—Ephemeris Time (ET). The 
pattern will be to discuss, in turn, the meaning, BpLODy, 
measurement and uses of ET; and since the “‘meaning’ 
of time is more properly a subject for philosophy, 
this area will be touched but lightly. Stress will be 
placed primarily upon the history, most recent methods 
of accurate time determination and some uses to 
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which this increased accuracy may be profitabl 
directed. 


2. Meaning and Definition of Ephemeris Time 


As to the origin of the name ‘‘Ephemeris Time, 
ET is derived from it being the independent variabl 
in the ephemerides, or time schedules, of celestia 
bodies; it is perfectly uniform in passage and is ofte 
used interchangeably with the name Newtonian time 
In Newtonian mechanics, time is looked upon a 
being instantaneously the same through the whole o 
space—a “now” which is present everywhere but i 
inexorably passing into a new. 

This view of time is to be contrasted with that a 
presented by relativity theory. Here, there is n 
longer the dichotomy of space and time, but rathe 
the concept of a space-time continuum. Time is n 
longer unconditionally the same for two observe1 
through all space, but is dependent not only upon th 
relative positions of the observers, but also. thei 
relative velocity. However, time will become Ney 


tonian for two observers if the relative velocity between 
hem becomes zero and will be essentially so for small 
alative velocities. 
‘It is enlightening to look upon Newton’s First Law 
. Motion (the Law of Inertia) as defining equal 
Bervals of time, i.e. uniform time. Paraphrasing the 
Law of Inertia, ‘Two intervals of time are equal by 
definition if a body subject to no forces moves over 
equal distances in them.” It is in this sense that we 
have so far used ET; however, let us look at a more 
practical, working neniuon of ET. 
_ At the September 1955 meeting of the International 
Astronomical Union (IAU) in Dublin, the IAU adopted 
HT as the measure of time to be used in precise work 
and, further, defined the unit of ET, the ephemeris 
second, as 1/31,556,925.975 of the tropical year 1900.0 
(1). ET is related to Universal Time (UT) by the 
Biuations HT = UT + AT and AT = 1.821 AL’y, 
where AL’, is the difference between the observed 
mean longitude of the Moon and the value tabulated 
in the new lunar ephemeris (Improved Lunar Ephemeris 
1952-1959). These equations will be considered as 
defining ET. 


3. Historical Developments Leading to the Rec- 
~ ognition of Ephemeris Time 

_ Before going further, it will be rewarding to look at 
the historical developments leading to so unique, 
singular a definition of the ephemeris second and to 
the relationship of ET and UT described. It will be 
further noted that there is a contending unit of time 
measurement? afforded by the atomic clocks developed 
m recent years; and it will be seen that these two 
basically different units of time? and the methods of 
determining them combine to give man his most 
accurate means ever for time determination. 

Data accumulated over several centuries point out 
the departures of the Moon, Sun and planets from 
their positions predicted from gravitational theory. It 
was noted early that the Moon’s motion had a secular 
component of acceleration; from records of eclipse 
studies, Hatiey concluded in 1695 that the motion of 
the Moon was being accelerated (2). It was later 
shown that a term of approximately 10’’T? (T in 
senturies) had to be applied to the mean longitude of 
the Moon to account for the acceleration; however, 
ypproximately half of the above term has been shown 
(0 be attributable to perturbative influences of the 
ther planets on the Earth-Moon system. 

In 1927, pe Srrrer found, from a study of lunar and 
“olar eclipse records covering about two thousand 
years, that the Moon’s secular acceleration (not 
ecounted for by gravitational theory) in mean longi- 
ude amounted to +5.22’” + 0.30” per century per 


1The atomic time unit has been termed the ‘‘essen’’ by 
x. M. CLEMENCE (4). 

2 The ephemeris second is of an astronomical nature and 
vased upon gravitation, whereas the essen is of a physical 
lature and based upon atomic processes. 


century (3). It was further noted by Cowxrtt in 1905 
that there was also a secular acceleration of the Sun’s 
motion (2). pm Srrrer’s value for the Sun’s accelera- 
tion is + 1.80’ + 0.16’. These secular accelerations 
have come to be attributed to the slowing of the 
EKarth’s rotation due to tidal friction. 

In addition to the secular variation mentioned, 
Newco, in 1870, found that there were also irregu- 
lar, possibly long period, fluctuations in the Moon’s 
motion which could not be accounted for by the 
current tables of the Moon (HANsEN’s Tables: pub- 
lished 1857) (2). Newcoms accounted for the fluctua- 
tions through the introduction of a theoretically 
unsubstantiated, empirical term which has become 
known as the Great Empirical Term. 

Irregular fluctuations in the Sun’s and_ planets’ 
motions were not recognized until after Brown (5) 
developed his lunar theory and constructed his Tables 
of the Moon. In the construction of his Tables (1919), 
Brown incorporated the Great Empirical Term and 
the secular acceleration term due to planetary pertu- 
bations, but he did not include the secular term due to 
tidal friction. Through the use of the fluctuation of 
the Moon, B, Spencer JonrEs was able to show in 
1939 that the observed departure from prediction of 
the Sun from Newcomp’s Tables (6) is expressible in 
the form ALg = + 1.00” + 2.97"T + 1.23”’T? + 
0.0747B (7); here, T is measured in Julian centuries 
from Greenwich Mean Noon 1900.0 (January 0.5 UT), 
and the time of observation is in UT units. B is defined 
as the fluctuation in the Moon’s longitude and is the 
difference between the observed mean longitude of 
the Moon and the tabular value from Brown’s Tables. 
Similar expressions may be derived for the several 
planets, where the secular and irregular terms are 
proportional to the mean motions of the planets. 

In the expression above, the first two terms are for 
adjustment to the epoch 1900.0. The third term 
represents the secular acceleration due to tidal friction 
slowing the Earth; however, note that the value is 
found to be different from the 1.80” value of DE SITTER. 
The fourth term is the irregular acceleration which is 
due to irregular fluctuations in the Earth’s rotation 
rate. It may at first sight appear strange that data 
from an observation of the Moon (B) would appear in 
the AL expressions of the Sun and planets. However, 
it must be remembered that we cannot detect the 
Earth’s fluctuations in rotation by purely terrestrial 
data—we must look to the celestial bodies for a “‘reflec- 
tion’”’ of these fluctuations in rotation rate. The ‘‘reflec- 
tion’? from the Moon is greatest, easiest measured and 
thus is used as the yardstick for the irregular fluctua- 
tions in the Earth’s rotation rate. 

It will be remembered that when UT was first 
instituted, it was thought to be Newtonian in charac- 
er; hence, it was thought that comparison of observed 
positions as functions of UT, and tabular positions, as 
functions of a perfectly Newtonian time (ET), for 
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celestial bodies would give zero residuals—except for 
uncertainties in constants, imperfections in gravita- 
tional theory and observational errors. 

As we have noted, the residuals are non-zero and 
are composed of secular and irregular components. 
Since the components of the residuals for the several 
planets are proportional to their mean motions, it is 
seen that the source of the residuals is most logically 
and rightfully assigned to the inconstancy of UT or as 
has been said, the inconstancy of the Earth’s rotation 
rate—since UT is measured by the rotation rate of 
the Earth. Had Newcomp’s Tables been entered with 
the time of observation in ET units, ALs would have 
been zero; thus, the ALsg is attributable to the difference 
in ET and UT. Where K = 248349/1.00” and is a 
measure of the rate at which the Sun’s mean longitude 
increases, we see that ET — UT = AT = K Als, or 
AT = + 249349 + 7253165T + 29s949T? + 1.821B. 
As before, the third and fourth terms represent, 
respectively, the secular and irregular components, 
and the first two terms are for adjustment to the 
epoch 1900.0. Hence, for the year 1900.0 (T = 0) 
ET = UT, in epoch and duration. To explain the 
source of the definition of the ephemeris second given 
earlier, we first note that the unit of UT is the mean 
solar second and is 1/86,400 of a mean solar day: 
hence, the number of mean solar seconds in the tropical 
year 1900.0 is 86,400 times the tropical year 1900.0 in 
mean solar days and is 31,556,925.975 seconds. Since 
ET = UT for the year 1900.0, the ephemeris second is 
equal to the mean solar second for that time, i.e. the 
ephemeris second is the fraction 1/31,556,925.975 of the 
tropical* year 1900.0. 

We have seen that ET — UT = AT, or ET = UT 
+ AT. It remains to show that AT = 1.821 Alyy. 
SPENCER JONES has shown that the observed departure 
from prediction in mean longitude of the Moon from 
Brown’s Tables is expressible in the form ALy = + 
4.65” + 12.96”T + 5.22”T? + B — Empirical Term 
(7). Here again the time of observation is in UT and T 
is in Julian centuries reckoned from the epoch 1900.0. 
Note that pm Srrrmr’s 5.22” value is retained for the 
secular component. Where K’ = 1”/1s821 and is the 
rate at which the Moon’s mean longitude increases, the 
Moon will increase its mean longitude by K’ AT in the 
interval AT; hence, the Moon’s increase in mean 
longitude in time AT is ALyar) = K’ AT = + 13.37” 
+ 39.71’T + 16.44”’T? + B. 

Thus, when the time of observation is referred to BT, 
the correction to Brown’s Tables becomes ALy — 
Alwar) = — 8.72” — 26.75"T — 11.22’’T? — Empiri- 
cal Term. It is principally through the inclusion of this 
corrective term that Brown’s Tables of the Motion of 
the Moon becomes the Improved Lunar Ephemeris 
1952-1959 (ILE). The effect has been to enter into the 

’’The sidereal year was first proposed, but it was later re- 


placed by the tropical year which requires no knowledge of 
the constant of precession and is thus more accurately known. 
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lunar theory the same unit of Newtonian time as used 
for the theories of the Sun and planets. | 

If now, the time of an observation in ET units is 
taken into the ILE, there will be no difference in ob- 
served and tabular values of mean longitude—except, 
again, that due to uncertainties in constants, imperfec- 
tions in theory and observational error. However, if 
the time of observation in UT units is taken into the 
ILE, there will be found an error AL’y in the observed 
and tabulated mean longitude which will be due to the 
difference in ET and UT for the epoch-of the observa- 
tion. Hence, ET — UT = AT = (1/K’) AL’y = 1.821 
Aa 

We have briefly seen the necessity of the transition 
from UT to ET and the essence of the principle of the 
conversion of the one to the other; however, in going 
from UT to ET we have given up a relatively accessible, 
variable time for a relatively inaccessible, invariable 
one. CLEMENCE (4) has indicated that, within a year, 
the year can be determined with a probable error of 1 
part in 10!° in UT units by astronomical methods using 
the Photographic Zenith Tube (PZT); however, to 
obtain the same accuracy in ET units, five years of 
observation would be necessary. Probable errors of the 
order of 1 part in 10° in ET units is expected from a 
single year of observation. 

In recent years, however, there have been developed 
atomic clocks which equal, if not exceed, the accuracy 
of time determination by astronomical methods—but 
without the lengthy delay. In 1955, Essmn and PARRY 
(8) announced that they had achieved an accuracy of 1 
part in 10° in time and frequency? determination 
derived from atomic processes and employing magnetic 
resonance techniques. There is thus available a new 
standard of frequency for calibration of other fre- 
quencies which permits of the same accuracy in a few 
minutes as given previously by astronomical methods 
requiring a period of a year. 

In the light of the Essmn-Parry article, BuLLARD 
(9) suggested that the astronomical second (ephemeris 
second) be abandoned in favor of a ‘physical second’ 
(essen) for the most refined experiments. With reference 
to the BuLLARD article, the IAU defined the unit of ET 
as previously stated. 


4. The Determination of Ephemeris Time 


Further response to the BuULLARD article came from 
CLEMENCE (10). Among other things, he pointed out 
that the atomic standard clocks could not establish an 
epoch. Astronomical observations of celestial bodies 
will be made to establish an epoch,and their motions 
will define the primary standard of time. The atomie 
clocks will serve as secondary standards to mark the 
passage of time between epoch determination—just as 
the quartz crystal clocks (and earlier still, mechanical 
clocks) have done in the past. 


“ Frequency is here used in the sense of the ratio of time 
intervals. 


We have noted that there is a long period measure- 
ment of ET using astronomical methods and a short 
period determination using atomic clocks. Taking 
these methods in turn, we remember that ET = UT + 
AT. UT is determined from observation using the 
PZT and AT is determined from observation using the 

Markowitz dual-rate Moon position camera. 


4.1. Long-Term Determination of Ephemeris Time 


4.1.1. Determination of Universal Time (UT). The 
PZT was devised by Ross (11) primarily to determine 
‘variations in latitude. The first was installed for its 
stated purpose in 1915. In 1933, it was adapted for the 
determination of UT and eventually replaced the 
transit instruments for this purpose. The PZT is 
directed toward the zenith® so that only those stars 
passing on or near the zenith will be photographed. The 
upper element of the instrument is called the rotary 
and contains the lens, photographic plate holder and 
drive mechanism. Light passing through the lens 
continues on through a fixed vertical tube until it is 
reflected by a basin of mercury situated separate from 
and below the tube. The light is reflected from the 
“pool of mercury back along the tube to the rotary 
where it is brought to focus on the photographic plate. 
~ Lrrreci and WItiis (12) have described two possible 
“methods of using the PZT for time determination: one 
method permits the stars to trail and the other provides 
drive to the plate holder permitting point images to be 
formed. Lirrexyi (13) later found the latter method to 
be best. In somewhat more detail, this method is as 
follows: The plate drive mechanism moves the photo- 
graphic plate in a west-east direction sufficient to 
balance the Earth’s rotational motion. Before transit, 
three, fifteen-second exposures of the star are made. To 
permit this, the photographic plate is moved approxi- 
mately 0.01 inch in the north-south direction between 
exposures. About five seconds before transit, the rotary 
is rotated 180 degrees; and starting about five seconds 
after rotation, there are formed three more fifteen- 
second exposures. 

Associated with each of the six star images are 
approximately twenty-seven time records, or eighty-one 
pairs total. Any pair, one before and one after rotation, 
would give the time of transit in units of clock time; 
but in practice, several pairs (a dozen or so) are used 
to obtain a more accurate value. Necessary correction 
to clock time for obtaining UT is determined from 
ephemeris data of the observed star’s right ascension for 
the date and the right ascension of the mean Sun for 
the epoch of observation. During a night’s observing, 
some fifteen to twenty stars are photographed, giving 
a probable error of about three milliseconds (4) in the 
mean result of mean solar time (UT) so determined. 
This method, over a period of a year, will give the 
accuracy previously stated of 1 part in 10'°. 


5 It will be noted that, the PZT being a zenith instrument, 
the effect of atmospheric refraction will be minimal. 


4.1.2. Determination of AT. In determining ET, AT 
must also be known in addition to UT. To determine 
AT, observations of the Moon’s actual position must be 
compared with its position given in the ILE; however, 
its actual position at a well-defined epoch is difficult to 
obtain using conventional techniques. The Moon’s 
greater brilliance in relation to the comparison stars 
results in an irradiation effect which makes the limb 
difficult to determine, and the Moon’s motion in the 
background of stars during the length of exposure 
creates an elongated lunar image whose center is 
difficult to ascertain. In addition to these difficulties, 
there is the further lack of a well-defined epoch for the 
observation. 

To overcome these difficulties, Marxowrrz (14) 
developed the dual-rate Moon position camera; the 
camera went into use at the U.S. Naval Observatory in 
1952. To eliminate the irradiation effect, the Moon’s 
image is intercepted by a plane, plate filter 1.8 mm 
thick which has a transmission factor of 0.001. A 
second, yellow filter having a hole in its center,® is 
located behind the first and is situated on the plate 
holder mechanism; its purpose is to cut out the blue 
rays of the stars. 

Normal tracking on the stars is not accomplished by 
the telescope drive, but rather the telescope is sta- 
tionary and the camera’s plate holder is driven in right 
ascension to permit tracking on the stars. To eliminate 
the effect of an elongated lunar image, the lunar 
filter is tilted at a rate which maintains the Moon’s 
image fixed in the background of stars. When the 
filters are parallel, the lunar filter passes the Moon’s 
rays undeviated; and at the instant the filters are 
parallel, a set of electrical contacts close to mark the 
epoch of the observation—in UT units. 

Reduction of the plate (14, 15) is as follows: Co- 
ordinates of 30-40 points on the Moon’s bright limb 
and about 10 comparison stars are obtained with a two- 
screw measuring engine. Since only three positions on 
the limb would be sufficient to define the center of the 
Moon, the method of least squares must be used when 
so large a number as 30—40 are used. The result is, then, 
the Moon’s center of figure given in right ascension and 
declination for the epoch of the observation. 

At this time, it is to be noted that, in practice, AT is 
not determined from AT = 1.821 ALM; but rather, by 
the method outlined in the ILE (16). Here, the value 
of right ascension of the Moon’s position is taken into 
the ILE and the value of ET is read out for which the 
value corresponds. This, then, is the epoch of the 
observation in ET units. The difference ET—UT for 
the epoch of the observation is AT. AT can also be 
determined using the declination, but it is not expected, 
however, that the two values agree exactly. The final, 
accepted value is obtained from an averaging of the 


6 The hole in the center of the second filter is to permit 
the Moon’s rays from the first filter to strike the photographic 
plate without further diminution. 
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two. A more accurate value of AT is obtained if the 
coordinates of the center of figure are first converted 
into those for the center of mass before entering the 
tables. 

From one observing station, AT should be deter- 
mined, over one lunation, with a probable error of Osl 
(15), and from observations over a period of one year, 
should determine absolute frequency to 3 parts in 10’. 
About twenty dual-rate cameras were scheduled for 
operation throughout the world during the recent, 
IGY; four are planned to be maintained in continued 
operation afterward. It was with consideration of four 
in operation that the value, previously given, of 1 part 
in 10°, over a period of a year, was given. 


4.2. Short-Term Determination of E’phemeris Time 


We have seen how man has gone from a time determi- 
nation employing the Earth’s rotation as a measure to 
one employing, essentially, the Earth’s revolution 
about the Sun. In all times though, it has been neces- 
sary to use the most accurately available Harth-bound 
timepiece for extrapolation between astronomical time 
determinations; however, these have been considered 
inferior timepieces, until quite recently, to the natural 
periodicity afforded by the celestial bodies and are 
daily checked upon, weather permitting. 

It is not to be overlooked, however, that there are 
advancements in the construction of more accurate 
time pieces; from the pendulum devices through the 
quartz crystal oscillators, we have arrived at the family 
of atomic clocks. We will look briefly at the ammonia, 
cesium and maser type clocks. 

4.2.1. Ammonia and Cesium Clocks: Resonance Ab- 
sorption Techniques. Fundamental to the ammonia and 
cesium type clocks, however, is the quartz crystal clock. 
Quartz is a piezoelectric crystal, and when properly cut 
and placed in a suitable oscillator circuit becomes a 
source of relatively stable electric oscillations. A 
representative value of these oscillations may be, say, 
100,000 cycles per second (eps). These oscillations are 
then counted-down through several cascaded multi- 
vibrator stages to a nominal value of 1,000 eps or, 
possibly, 60 cps. Here they undergo a few stages of 
power amplification and eventually drive what may be 
an ordinary electric clock. The fundamental oscilla- 
tions could as well have undergone frequency multipli- 
cation and amplification and eventually been broadcast 
as radio signals. However, a drawback to the quartz 
clock is a secular change in frequency with age. 

The ammonia (17, 18) and cesium (8, 17) clocks are 
similar in that they provide a means for stabilizing the 
’ frequency of the quartz oscillator by comparison to an 
atomic or molecular standard of frequency. The method 
of comparison is one of resonance absorption. In the 
case of the ammonia molecule, the frequency of ab- 
sorption used for comparison arises from the nitrogen 
atom oscillating between its two possible positions 
with respect to the three hydrogen atoms. The three 
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hydrogen atoms form the base of a triangular pyrami 
whose apex is occupied by the nitrogen atom; however 
the nitrogen atom can oscillate “through” the base to 

second possible apex. 

For cesium, which is an alkali metal and hence 
hydrogen like in structure, the absorption frequency} 
arises much as does the 21 cm absorption of the neutra 
hydrogen atom. The spin of the outermost, valen 
electron can be either parallel or antiparallel to the 
nuclear spin, giving rise to two distinct energy states 
ot the cesium atom. In the case of the cesium atom, the 
frequency of absorption used for comparison arises 
from the difference of energy of these two states. * 

In both the ammonia and cesium clocks, quartz 
crystal oscillations are multiplied in frequency until 
they are equal to the resonance absorption frequency of 
the atomic standard. Then, in the ammonia clock, the 
signal is applied to one end of a waveguide filled with 
ammonia gas. At the farther end of the waveguide is a 
microwave detector. If the input to the detector is not 
a minimum (absorption in the waveguide a maximum), 
the detector output will alter the quartz crystal oscilla- 
tor frequency until it is a minimum. 

In the cesium clock, the method is somewhat differ- 
ent. Here, the multiplied oscillator output reacts with a 
beam of cesium atoms to effect absorption of the 
microwave energy; the effect of the absorbed energy on 
the beam is to increase the number of the cesium atoms 
in the higher of the two energy states discussed. After 
reaction with the microwave energy, the beam of 
cesium atoms passes through a magnetic field which 
isolates the atoms in the different energy states and 
permits those atoms in the higher energy state te 
continue on to a detector. If the detector output current 
is not maximal, the crystal oscillator frequency will be 
shifted until it is. 

4.2.2. Maser Clock. The maser (18, 19) differs from 
the ammonia and cesium clocks primarily in that it 
does not depend upon the mechanism of absorption. 
Instead, the frequency standard of the atom itself is 
used directly—much as the initial crystal oscillations of 
the quartz clock had been. The maser can use any of 
several different atomic or molecular sources; but, fot 
illustration, we will consider only the use of the cesium 
atom. 

In a group of cesium atoms, there will be some atoms 
in each energy state previously discussed. Forming the 
atoms into a beam and passing them through a mag: 
netic field (as was done in the cesium clock), a bear 
can be generated which is predominantly made up ot 
atoms in the higher of the two energy states. If such ¢ 
beam, of sufficient intensity, is fed into a resonant 
cavity of the proper dimensions,’ the cesium atoms wil 
return to the lower energy state and give up thei 
energy to the cavity. The cavity oscillations are ther 


‘The dimensions of the cavity are such that it is resonan’ 
to the frequency of oscillation given off by the atoms whe! 
returning to the lower energy state. 


upled out by means of a waveguide to the ampli- 
ying and count-down circuits. If a beam of intensity 
nsufficient to sustain oscillations in the cavity is 
resent, an input signal frequency equal to the cavity 


requency will stimulate the cesium atoms to radiate, 
with the result that the input signal frequency appears 
increased, or amplified. This is the principle of the 
maser as an amplifier. 


_ The accuracy of the ammonia and cesium clocks are 
Bendent upon the sharpness of the associated absorp- 
tion curve; the cesium clock overcomes the Doppler 
b1 Baching which is inherent in the ammonia clock and 
is thus the more accurate of the two. It was using a 
eesium clock that Essen and Parry (8) achieved their 
spoch-making frequency determination. The maser 
presently holds the highest hopes for the greatest 
necuracy; it is anticipated that the maser may achieve 
an accuracy of 1 part in 10% when using the cesium 
atom (17). 


5. Some Uses of Ephemeris Time 


_ When speaking of the ‘‘uses” of an improved unit of 
ime, one is faced with an enormous problem of selec- 
ion; for in all scientific and technical areas of human 
endeavor, a more precise, uniform standard of time 
determination will find immediate use—if for only 
pushing the value of some constant out to another 
significant figure. But it is just this knowledge of the 
added figure that may open a new corridor into the 
mknown. 

In astronomy, an effect of the use of the new time 
ET) has been the removal of a glaring inconsistency 
tween lunar theory and planetary theory—they now 
se the same time variable in the ephemerides. If now, 
inaccounted-for errors exist between observation and 
heory for the position of the Moon, Sun and planets, 
here leads a way to improved gravitational theory. 
Discrepancies may indicate necessary revision to orbit 
barameters or masses; in general, the new time should 
bermit improvement in both lunar and _ planetary 
heory. 

The physicist and radio engineer will use the new 
ime unit for expressing frequency much as previously 
put to a higher accuracy. Where fy is the frequency 
based on UT units, fp is the frequency based on ET 
nits and D AT is the change per annum of AT, fy = 
gi = AT/31,556,925.975) = fy(i — 3.1689 X 
O0-DAT) (1). fy is determined by comparison to the 
jational time services radio broadcasts and the D AT is 
ublished regularly by the time services. Once the fx of 
he atomic frequency standard (e.g., the absorption 
requency of the ammonia and cesium clocks) has been 


8 Maser is a shortened expression for microwave amplifica- 
ion by stimulated emission of radiation. 


sufficiently accurately determined, ET may be broad- 
cast directly by the time services. 

This (the determination of the essen in ET units) 
leads to one of the most fundamental questions which 
may be answered by the new time measures: Is the 
ephemeris second and the essen measuring the same, or, 
is the time from gravitational measurements the same 
as that from atomic measurements? In kinematic 
relativity (4, 9), there are two natural time scales 
present, which do not keep in step. Minne, Dirac and 
others have suggested they be associated with atomic 
and astronomic frequencies. A much looked-for result 
is the manner in which the essen and ephemeris second 
keep in step. As has been often and aptly said, ‘Time 
will tell’—the answer to this question and many more. 
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Abstract 


Calculations are made, for the particle spectra in the most 
intense regions of the radiation belts, of the attenuation of the 
particle flux by shielding on a space vehicle, and of the bio- 
logical dose rate due to the flux which has penetrated the 
shielding. Both protons and electrons are considered in the 
inner belt, but only electrons in the outer belt. Secondary 
particles and bremsstrahlung produced by these particles on 
interacting with a shield are also considered. It is shown that, 
despite the much higher electron fluxes, the protons are of 
much greater biological significance. 


Introduction 


With the discovery of trapped radiation belts around 
the earth at altitudes between approximately 1,006 
km to 50,000 km, there has arisen the question as to 
the magnitude of radiation dose that would be received 
by the occupant of a space vehicle which passed through 
these regions. In addition, there has been raised the 
question of the amount of shielding that would be 
required in a space vehicle to protect the occupant, 
should the dose rates in the radiation belts prove to be 
high or the time spent therein long. 

It is the purpose of this paper to discuss quantita- 
tively the question of dose rate and effectiveness of 
shielding in reducing dose rates in the most intense 
regions of the radiation belts. Estimates are made, 
based on experimental data, of the particle spectra in 
these regions. Calculations are then made of the at- 
tenuation of incident particle flux due to various 
amounts of shielding, and of the biological dose rate 
due to the flux which has penetrated the shielding. 

The next four sections are concerned with the 
inner belt protons, which are known to be the most 
penetrating particles in the radiation belts, and with 
the secondary radiation generated by the protons on 
interacting with a shield. The last two sections are 
devoted to the electrons in the inner and outer belts, 
respectively, and to the bremsstrahlung which they 
produce. 


Incident Proton Spectrum 
With data now available [1, 2, 3, 4] regarding the 
flux and energy distribution of the protons of the 
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inner radiation belt, it is possible to make reasonabl 
estimates of the incident radiation intensity above 4 
Mey in this region around the earth. Below 40 Me 
theoretical forms of the energy spectrum [5] exist 
based on neutron leakage from the atmosphere. Th 
total omnidirectional flux above 40 Mev has bee! 
given by Van Allen [2, 4] as 2 X 10* protons/cm?-se 
in the heart of the inner belt. Some doubt [6] exists a 
to the validity of the Freden and White spectrum [8 
over the whole belt region, but it certainly represent 
an average over that part of the belt through whiel 
their emulsions passed, and no better form for thi 
spectrum is yet known. 

Figure 1 shows differential energy distributions 0 
inner belt protons which fit the experimental dati 
now known. P; is the form given by Hess [5] belov 
100 Mev, joined to the form exp (—#/120) abov 
100 Mev [3], having a total of 2 X 10* protons/cm?-se 
above 40 Mev. P, is an alternate spectrum whicl 
fits the data: 


P(E dE; 
= 1.04 X 10!E;' dE,(40 < E; < 100 Mev) a 
P(E dE; ; 
= 1.04 X 10°H;° dE\(E; > 100 Mev) 


The subscript 7 refers to incident particles, and capita 
P refers to the absolute omnidirectional flux distribu 
tion. Lower case p is used to designate a spectrun 
normalized to one incident proton/em?2-sec total flw 
of energy above 40 Mey. po(#;) is plotted in Figur 
2, along with spectra calculated later. By using Vat 
Allen’s total flux, PZ) = 2 X 10*p(#). The biologiea 
dose calculations reported here are based on the forn 
po; )dE; for the incident spectrum. 


Proton Attenuation and Spectrum Degradation 


The incident protons impinge upon the outer shel 
of a vehicle, and are partly absorbed and degraded t 
a new spectrum by it. It is then assumed that the 
enter some absorber designed to reduce the dose insid 
effectively. Polyethylene (CH2) has here been chose 
as the material for this absorber because it effectivel: 
stops protons, on a weight basis, as well as or bette 
than most other substances, with a minimum ¢@ 
nuclear reaction products. The radiation entering th 
vehicle interior, i.e. the astronaut’s compartment 
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_ Fra. 1. Inner radiation belt proton differential energy spec- 
tra. Circles are the experimental data of Freden and White 
[3]. P: is one of their curves fitting the data. P: is the power 
law spectrum used in this paper, Equation (1). 


thus consists of degraded energy incident protons, 
neutrons and gamma rays generated in the metal 
shell by (p, n) and (p, y) reactions, plus penetrating 
electrons and electron-produced bremsstrahlung. These 
will be discussed in turn. 

Tables of range vs. energy for protons in various 
materials [7] were used to calculate the energy de- 
eradation as the protons pass through matter. Using 
only R.(/), range as function of energy in a material 
Z, one may write: 


RAE;) = AG 45 R,(E;) (2) 


where X. is the thickness of the absorber Z, and HL; 
is the final energy of the proton after passing through 
this thickness. /; is the incident proton energy. 
The form of R.(/) may be approximated by a 
power law, 
RAH) = KE", (3) 


where both & and n depend on FE as well as (slightly) 
upon Z. If we denote initial and final energies by sub- 
scripts 7 and f, and consider only a single material, 
then Equation (1) becomes: 


i ae ea ae (4) 
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Fie. 2. Proton differential flux spectra, normalized to unit 
incident flux over 40 Mev. The top curve is the unit incident 
flux, p. (#). The other curves are degraded energy spectra 
resulting from passage of the incident flux through the indi- 
cated shielding. 


Differentiating, 
pis ky N; E; es 


kyng Eyl 


dH; i; (5) 
Hence, if p(H;)dE; particles enter the absorber 
having energies in the range H; to H; + dH;, where 
FE is large enough so that traversal takes place through 
a thickness X, then this same number will emerge 
having energies in the range H; to H + dH;,, neglecting 
straggling. In equation form, 
p(E;) dk; = p(#;) dk; 


kms Ey?” WE (6) 
ke ee 
Over the range of interest here, 40 to 1000 Mev, and 
for relatively small energy losses in the absorbers, k 
and n may be considered constants, and the ratio 
kini/kyny set equal to one. 

In (6)-#; is related to H; in such a way that (2) 
is satisfied if X is thin enough so that #; is not entirely 
absorbed. This is equivalent to integrating the energy 
loss rate equation of Bethe-Bloch [7] over the absorber 
thickness and subtracting it from the original energy. 

The recipe for obtaining LH; from H; and X, is as 
follows: for H;, get R.(H;) from the range-energy 
table, and subtract X. from it. Reenter the table to 
find Hy corresponding to this difference, R,(H;). In 
the conversion of p(H;) to p(H;), the value of n for 
polyethylene was taken as 1.81 up to several hundred 
Mev,* and for steel, n = 1.77. The calculation of the 


= p(H;) 


* Actually n decreases slightly with increasing H#, and k 
increases. The error in p(#) caused by neglect of these varia- 
tions is slight, and results in a few percent overestimate in 
dose. 
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final spectrum p,(H,) resulting from passage of protons 
of original spectrum p.(H#;) through various amounts 
of steel and polyethylene is a matter of setting up 
tables to follow this recipe (2) in conjunction with 
(6). Degraded spectra resulting from such tabulated 
calculations are plotted in Figure 2, along with the 
incident spectrum. 

The degraded spectra of Figure 2 have been cal- 
culated assuming the protons penetrate the shield 
normally. In actuality non-normal penetration will 
also take place. Because of this slant penetration, the 
effective shield thickness will be greater than the 
actual value, resulting in somewhat lower (perhaps 
10 to 30%) fluxes than calculated here. 


Nuclear Reactions 

As protons pass through the steel shell, or other 
metal structure on and in a space vehicle, some nu- 
clear interactions will take place. These will be (p, mn), 
(p, y) and (p, p’y) events, producing secondaries which 
are capable of delivering radiation doses to the as- 
tronaut. Hence some consideration was given to these 
effects. Data from various sources were used to estimate 
cross sections for these processes as a function of 
energy from a few Mev up to 1000 Mey. Both iron 
and molybdenum were considered. 

The data of Crandall and Millburn [8] were used, 
along with inelastic scattering cross section data of 
Hicks and Kirschbaum [9], and others [10, 11, 12], 
to arrive at the following conclusions regarding neutron 
production at incident proton energies over 40 Mev: 
(1) in steel, 0.07 X neutrons per incident proton, and 
(2) in molybdenum, 0.15 X neutrons per incident 
proton, where X is the thickness of the vehicle shell in 
em, and the protons are distributed in energy ac- 
cording to Equation (1). Consideration of energies 
below 40 Mev leads to neutron numbers several orders 
of magnitude less than the higher energy case. 

Assuming an approximately constant cross section 
of 1 barn (=10~*4 em?) for inelastic collision excitation 
of nuclei in the vehicle shell, the number of gamma 
rays generated is 4 X 10°? X photons/em? (of a few 
Mev from many nuclear levels) per incident proton /em?2 
where again X is the shell thickness in centimeters. 

These neutron and photon energy fluxes are much 
smaller than those due to direct penetrating protons 
for most feasible shielding configurations. 


Dose Rate Calculations 

The proton flux p(Z) calculated earlier and plotted 
in Figure 2 for a given shield represents the differential 
energy distribution of normally-penetrating protons 
inside an empty vehicle. This is the flux that would be 
measured by a radiation detector small enough that 
self-shielding effects are negligible. The absorbed dose 
rate in this small detector is related to the differential 

particle flux p(#) dH by 


; dE 
= = 576 X 10° ‘f ome 


- (EH) dE rad/hr (7) 
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where dH/dx (E) is the energy loss rate in the detect 
in Mev-cm2/gram, and one rad is 100 ergs/gram. 

Equation (7) may be used to calculate the ener 
absorption rate at any point in a mass of materi 
provided that p(#) is known at that particular pot 
and dE/dx as a function of energy is known for th 
material. Hine and Brownell [13] give data from whi 
one may obtain dE/dx (E) for tissue. This tissue energy 
absorption rate is plotted as the lower line in Figure 3 
On using this curve in Equation (7), along with thi 
appropriate p(#) curve from Figure 2; one obtains 4 
absorbed energy dose rate in a small detector as_ 
function of vehicle shielding thickness. These vesuill 
in tissue rad/hr per unit incident flux, are tabulated i 
the second column of Table I. 

The effectiveness of radiation in causing biologian 
damage depends not only upon the amount of radian 
energy absorbed, but also upon the absorbed energy 
density and upon the biological effect (‘end point” 
being considered. The absorbed energy density 
terms of linear energy transfer, kev/micron, will bi 
large for slow-moving heavy particles such as proton 
of a few Mey. Consequently their effectiveness 1 
enhanced with respect to particles with lower specifi 
ionization. 
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TABLE I 
Inner Radiation Belt Dose Rates (Protons Only) 
Small Detector Proton Dose Rates 

Shielding: Steel Body Body 
 - CHy Surface Center 
a 3 Woitineident Aas For 2 X 104 Dose Rate|Dose Rate 
4 p/cm? sec 
— 
e (Inches) afl sani rad/hr | rem/hr rem/hr rem/hr 
0 + 0 3.5 22 70 440 220 a0 
0 se 0.47 3.0 9.3 60 30 6.7 
Mees 2 © 9) °0.29)) 1.7 |.5.8 | 34 17 5.8 
0 + 4 0.16 0.83 | 3.2 16.6 8.3 4.5 
0 + 8 0.09 0.38 ist 7.6 3.8 3.5 
oo.04 + 0 1.25 8.5 25 170 85 6.6 
0.04 + 0.5 | 0.36 2.3 22 45 22.5 6.1 
~0.04 + 1 0.25 aS 4.9 30 15 5.7 
0.04 + 2 0.19 1.0 Dretl 20 10 4.8 
mp.04 -+ 4 0.11 0.50 | 2.1 10.4 5.2 3.6 
0.04 + 8 0.06 0.30 igi 4.8 2.4 2.2 
10: -- 0 0.5 2.5 10 50 25 6.7 


Relative Biological Effectiveness (RBE) of a radia- 
tion type is defined as the ratio of the absorbed dose of 
a standard x- or y-ray quality to the dose of that 
radiation type required to elicit the same biological 
effect. A full discussion of the RBE concept and the 
difficulties with which RBE determinations are fraught 
is out of place here [14, 15]. 
~ In establishing maximum permissible doses for man, 
and in safety considerations generally, the RBE for 
heavily ionizing particles is taken to be 10 or 20. For 
protons up to about 10 Mev, an RBE value of 10 
is given by Hine and Brownell, and this value has been 
taken for the calculations here. In the range 10 Mev 
to 100 Mey, RBE for protons is here assumed to 
decrease gradually to 1.0, and from 100 Mev upwards, 
it is assumed to remain constant at 1.0. The upper 
curve of Figure 3 shows the assumed rem(£), that is 
RBE  X rad, per unit incident flux, as a function of the 
proton energy. In a military situation where safety 
considerations might give place to strategic ones, an 
RBE value of 1.0 over the whole energy range might 
be used, that is, the lower curve of Figure 3. This is 
equivalent to increasing the allowable dose by a factor 
of about 5 or 6. 

The biological dose rates corresponding to the 
energy-absorbed dose rates of Table I, Column 2, have 
been computed by evaluating the expression: 


oe = 3600 [ p(£) rem (EH) dE rem/hr. (8) 
0 


These dose rates, in rem/hr per unit incident flux, are 
tabulated in Column 3 of Table I. The ratios of the 
yalues in Column 3 to the corresponding values in 
Column 2, give the effective, or average, RBE value 
for the particular shielding thickness in question. The 
effective RBE decreases as shielding thickness in- 
creases. This is consistent with the shift of the peak 
of p(#), Figure 2, to higher energy as thickness in- 
sreases. Columns 2 and 3 of Table I are applicable to 


unit incident flux (above 40 Mev), as defined earlier. 
Columns 4 and 5 give the corresponding total proton 
dose rates to be expected in the heart of the inner 
radiation belt, based on present data. Some of these 
calculated data are also presented in Figures 4 and 5. 

The preceding calculations in this section have been 
for the case of a small detector, in which self-shielding 
or self-absorption effects are negligible. For an as- 
tronaut in a space vehicle, self-shielding of one part of 
the body by other parts is of considerable importance, 
in that the dose rate at any given point in the body will 
be less than the dose rate measured by a small detector 
at the same location inside the empty shielded vehicle. 
The exact calculation of the dose rate as a function of 
depth in tissue would require a detailed knowledge of 
the orientation of the body and its location inside the 
vehicle shield, and would involve computing a p(/) 
distribution at each point in the body. Such a detailed 
calculation has not been made here. Instead, to make 
the calculations tractable, we have assumed that the 
human body can be dosimetrically represented by a 
sphere of tissue twelve inches in diameter, and that 
this sphere is located near the center of the shielded 
volume. 


SMALL DETECTOR DOSE RATES — INNER BELT 


POLYETHYLENE ONLY 


.04” STEEL PLUS 
POLYETHYLENE 


DOSE RATE RAD/HR 


PROTONS ONLY 


ELECTRON 
BREMSSTRAHLUNG 


POLYETHYLENE ONLY ONLY 


“0 2 4 6 8 10 12 14 
INCHES POLYETHYLENE SHIELDING 

Fra. 4. Absorbed dose rates in rad/hr due to protons and 
electrons in heart of inner belt, for various polyethylene thick- 
nesses with and without 0.04” steel vehicle shell. For the region 
0 to 2” for the case of 0.04” steel and polyethylene, and to 1” 
for polyethylene only, the shape of the dose rate curves has 
been estimated as indicated by the broken lines. 
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BIOLOGICAL DOSE RATES 
INNER BELT 


DOSE RATE IN REM/HR 


BODY SURFACE 
DOSE RATE 


BODY 
CENTER 
DOSE RATE 


—-—-—-— PROTON DOSE RATE 


TOTAL DOSE RATE 


0 2 4 6 8 10 12 
INCHES OF POLYETHYLENE SHIELDING 


Fic. 5. Biological dose rates, in rem/hr, due to protons and 
electrons in heart of inner belt, for various polyethylene thick- 
nesses, with and without 0.04” steel vehicle shell. Dotted lines 
are dose rates due to protons only. Solid lines are total dose 
rates due to both electron bremsstrahlung and protons. For 


a1” 


the region 0 to }” for the case of 0.04” steel and polyethylene, 
and to 1” for polyethylene only, the shape of the body center 
dose rate curves has been estimated as indicated by the broken 
lines. 


Considering the omnidirectional character of the 
incident flux, and the proton absorption in the body 
which gives self-shielding, we estimate that the dose 
rate on the surface of this tissue sphere is one half the 
dose rate in a small detector placed at the same point 
with the sphere removed. These estimated body surface 
dose rates are tabulated in Column 6 of Table I. 

It is possible to calculate the dose rate at the center of 
this idealized spherical body. In order to use Equation 
(8), p(#) must be known at the point where the dose 
rate is desired. This distribution was obtained from 
those plotted in Figure 2 by further degrading them 
through six inches of tissue. lor each vehicle shielding 
configuration, the resulting p(/) at the center of the 
sphere was folded into the rem(#) curve and nu- 
merically integrated. The biological dose rates thus 
valculated are listed in Column 7 of Table I after 
multiplying by the factor of 2 « 10‘. These figures are 
correct for the geometry and other conditions assumed, 
and reasonably close for the actual situation, the 
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human trunk. These dose rates are also shown 
Figure 5. 

A more detailed computation has begun, using 
general incident proton flux parameters, and for various 
other shielding materials and configurations. As the 
spectral form of the proton flux becomes better kno W 
more accurate dose calculations will thus be possible 


Electron Doses (Bremsstrahlung) 


So far only protons have been considered. In the 
inner belt, omnidirectional electron flux-values of 10 
em~ sec-! with energy above 20 kev, and 10’ cm™ 
sec! with energy above 600 kev have been given [2, 4]. 
These data are consistent with a differential spectrum 


electrons (9) 


9 —7.95E 
N.(#) dE = 93 X10" "dE 


which is used here for want of more definitive spectral 
information. Less than 1% of the electron energy flux 
can penetrate a steel shell as thick as 40 mils. For the 
thick shell case, the bremsstrahlung energy flux pro- 
duced in the shell by electrons in the energy range 
Eto E + dE is: 


ZE 
— 2 Wh ae 10 
top a N.(E) (10) 
or, for steel, and the flux (9), 
r dE = 3 X 10°E* exp (—7.95E) de “_. GB 
cm*-sec 


The total bremsstrahlung energy is thus: 
[ lr dE = 12 X 10° Mev/em™sec. (12) 
0 


Dividing by 10° incident electrons/cem?-sec gives a 
bremsstrahlung energy flux per ‘‘unit incident electron 
flux” of 1.2 X 10-* Mev/cm?-sec. 

The preceding calculations have made two assump- 
tions concerning bremsstrahlung generation, both of 
which result in some overestimate of the x-ray doses. 

1. The ratio of bremsstrahlung production to ioniza- 
tion loss by electrons, and approximately to total 
electron energy loss, is given by HZ/800. This ratic 
actually apples to rates of energy loss, but is here 
used for total energy loss. 

2. The direction of the incident electron flux is the 
direction of the resulting photon flux. 

As this photon flux proceeds into the vehicle, it wil 
be attenuated by any further shielding used. The 
maximum of the distribution function (11) occurs at 
0.25 Mev, so a choice of this energy was made as the 
most probable photon energy. To calculate the reductior 
in the astronaut’s absorbed dose effected by shielding 
inside the vehicle shell, we have used ‘poor geometry’ 
attenuation, since scattered photons may also con: 
tribute to dose. Alternatively, the total (“good geom. 
etry”) attenuation coefficients could have been used 
along with build-up factors. The poor geometry attenua 
tion coefficient for CH» at 0.25 Mev is 0.03 cem?/gram 


These considerations lead to dose attenuation factors 
of 0.93, 0.87, 0.76, 0.58 respectively for 1”, 2”, 4”, and 
8” of Pole thelene Applied to the ponerated x-ray 
energy flux of 1.2 X 10° Mev/cm?2-sec or 2.4 roentgen/ 
hour, the astronaut’s absorbed tissue dose rate is 2.1 
vad /hr. The small detector dose rates due to brems- 
‘strahlung are given in Table IJ, Column 2, and in Figure 
4. The body surface dose rate is estimated as three- 
fourths of the small detector rate, because of self- 
shielding in the body, and is given in Column 3. The 
body center dose rates, calculated by considering the 
attenuation of the x-ray flux in six inches of tissue, 
are listed in Column 3 of Table II. The contribution of 
the inner belt electron bremsstrahlung to the astronaut’s 
dose is shown in Figure 5, and is the difference between 
the solid and dotted curves. 


Outer Belt Electrons 


_ Estimates of the outer belt electron fluxes have 
also been made [2, 4]. These electrons are generally 
lower energy, there being 10" electrons/cm?2-sec, having 
# > 0.02 Mev, and 10° electrons/cm?-sec, having 
E > 0.2 Mev. The number of electrons of energy 
greater than 4 Mev is negligible. The data given are 
consistent with the differential distribution: 


128/028 oy ay electrons 


eet 13) 
cm*-sec 


N.(E) dE_= 83 X 10°e 


: Using the method outlined in the preceding section, 
by substituting (13) into (11), the total bremsstrahlung 
energy generated in a steel shell of thickness greater 
than a few mils is found to be 9.5 X 10° Mev/cm?-sec. 


The electron energy for which the maximum brems-. 


strahlung energy flux occurs is 50 kev, and so this 
energy was chosen as the most probable photon energy 
for purposes of calculating absorbed dose and attenua- 
tion. 

If one assumes that this 50 kev average energy 
bremsstrahlung flux is generated near the outer surface 
of the 0.04” steel shell, then the attenuation and filter- 
ing effect of the steel must be considered. Using 1.1 
em?/g as the ‘“‘poor geometry” attenuation coefficient 
of steel 0.04” will reduce the flux by a factor of two. 
This underestimates the reduction of dose because of 
the filtering effect of the steel which strongly absorbs 
the lower energy x-rays by photoelectric processes. 


TABLE II 
Inner Belt Electron Bremsstrahlung Dose Rates 


Shield: Small Detector Body Surface Body Center 


Steel + CH2 Dose Rate Dose Rate Dose Rate 
(inches) ae ee elec/ rad/hr (= rem/hr) | rad/hr (= rem/hr) 

0.04 + 0 2.1 1.6 1.4 

0.04 + 1 1.9 1.5 1.3 

0.04 + 2 1.8 1.4 1.2 

0.04 + 4 1.6 iL. Il 

0.04 + 8 12 0.9 0.7 


For no steel (CH: only) multiply above figures by 6/26. 


At 50 kev, the poor geometry absorption coefficient 
for CHp is 0.03 em?/gram, and that for tissue was taken 
as 0.04 cm?/gram. Hence, if no interior shielding is 


provided, the entire bremsstrahlung flux emerging 


from the steel shell, namely 4.8 X 10° Mev/cm2-sec 
will impinge upon the astronaut’s body. The resultant 
“small detector dose rate” is 11.5 tissue rad/hour. The 
body surface dose rate would be approximately three 
fourths of this amount. The center dose rate in the 
idealized spherical body, that is the dose delivered by 
this flux after passing through six inches of tissue, is 
6 rad/hr. 

If polyethylene is used for interior shielding, the 50 
kev x-ray dose attenuation factors for 1”, 2”, 4”, and 
8” are, respectively, 0.91, 0.82, 0.67 and 0.45. These 
factors are to be applied to each of the dose rates 
mentioned, i.e. small detector, body surface, and body 
center. 

The dose rates given here are conservative from the 
point of view of safety. Lower dose rates would result 
from the choice of lower atomic number materials for 
the outer shell. For very thin shells, and no interior 
shielding, some electrons might penetrate, which could 
cause skin burn doses. However, in view of the inner 
belt, shielding requirements, it is likely that these 
direct effects will be eliminated. 

Because of current doubts concerning the outer belt 
electron total flux, it is best to consider these dose 
rates in terms of a unit incident flux having the spectral 
form (13). This is obtained by dividing the figures 
given by 10". Thus the small detector dose rate is 
1.2 X 10~, the body surface dose rate is 9 X 107", 
and the body center dose rate is 6 & 10-" rad/hr for 
unit incident flux. 

Lack of detailed information regarding the outer 
belt electron flux energy distribution has justified the 
types of approximation made here. Refined calculations 
to eliminate the bremsstrahlung generation and ab- 
sorption approximations are being made so that, 
when better flux information becomes available, it may = 
be introduced into these computations to obtain- 
more exact dose rates. 


Conclusions 


We have calculated the biological dose rates inside a 
shielded space vehicle due te the protons and electrons 
which are present in the most intense regions of the 
radiation belts. The protons in the inner belt are by far 
the largest contributor to the dose, due to their high 
energies and relatively great biological effectiveness. 
The contribution to the dose of secondary radiations 
due to protons interacting with the shield is small, 
particularly for thinner shields: The electron-produced 
contribution is also relatively-small, despite the large 
incident electron flux and the simplified bremsstrahlung- 
production analysis employed here, which overestimates 
the bremsstrahlung to be expected. 

Another noteworthy point is the rapid decrease of 
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dose rate, from the unshielded case, with the addition 
of small amounts of shielding. For thicker shields a 
given increment of material is less effective in reducing 
dose rate. This effect can be seen from the curves in 
Figures 4 and 5. From an engineering and cost stand- 
point, it would appear better to provide a small amount 
of shielding and to select vehicle trajectories so that 
not too much time is spent in the heart of the inner 
radiation belt. 

While the dose rates are lower for the thicker shields, 
this absorbed energy is distributed more uniformly in 
the astronaut’s body. For the thinner shields, a large 
fraction of the total dose is absorbed at the body surface. 
Thus a long-term exposure in the inner radiation belt, 
which would require a thick shield, will result in near- 
uniform absorbed dose distribution in the body. A 
short-term exposure, as might result from rapid passage 
through the belt, and which would permit thinner 
shields, would give a larger fraction of the absorbed 
dose to the body surface. These are additional facts 
to consider in space trajectory planning. 

A more detailed calculation is currently being pro- 
grammed for automatic computation of biological 
doses. General parameters specifying radiation belt 
fluxes and spectra as function of altitude, and general 
shielding parameters specifying atomic numbers, den- 
sities, and thicknesses will be used. Bremsstrahlung 
calculations will be refined. A body depth-dose calcula- 
tion will be made. This program should lead to more 
exact dose rates based on the more reliable particle 
flux information which is sure to be available in the 
future. 
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The Influence of Heart Action and the Cir: 
culation of the Blood on Manned Satellite 
Attitude Control’ 


Lawrence R. Zeitlin? 


Introduction 

Among the internal disturbing forces which must be 
considered in the design of an attitude control system 
for a manned satellite are those caused by the physio- 
logical processes of the operation, particularly those 
produced by heart action and the circulation of the 
blood. 

' This work was sponsored by the Flight Control Laboratory 
through a subcontract under Contract AF33(616)-6674 with 
Systems Corporation of America. 

? Dunlap and Associates, Ine., Santa Monica, Calif. 
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l. Ballistic 
other impulsive forces. 

2. Mass displacement—caused by pooling of blooc 
during sleep, digestion, ete. 


3. Torques or couples—caused by dissymmetry 0 
blood flow. 


heartbeat anc 


Ballistic or Recoil Forces 


It has long been recognized that as a gun kicks wher 
it 1s fired, so also does the body recoil in response to th 


~~ 
K 


3 Fia. 1 


t ropelling and impact forces which move the blood in 
a vascular system. 

5. The recoil motions of the body were first described 
by Gordon in 1877 and elaborated upon by Henderson 
in 1905, the honorary godparents of the modern 
science of ballistocardiography, the study of heart 
function by the measurement of gross body movement 
(5). 

Essentially, a ballistocardiograph is a spring-sus- 
pended bed upon which the subject is tightly strapped. 
Measuring devices pick off longitudinal motion and 
amplify it for graphic presentation. When plotted 
against time, the motions of the ballistocardiograph 
produce a graph of the general form shown in Fig. 1. 

H represents auricular ejections; J represents ventric- 
ular systole, headward ejection; J represents footward 
motion of blood caused by turn of aorta; and K repre- 
sents the slowing of the downward moving blood. 

The equations describing this motion are derived 
from the principle of conservation of momentum, 1.e., 
mass of ejected blood X velocity of blood = (mass of 
body + mass of ballistocardiograph bed) X velocity 
(see Fig. 2). 

For a standard low frequency ballistocardiograph 
the movement (D) of the ballistic system in em. be- 
tween the J and J peaks can be expressed as: 


980F 
An? foe (Mi + Mz) 


D= 


where fp = natural period of the ballistic system 
F = force in grams required to move system (1, 
fGR 
If we consider the vehicle and its occupant as an 
extremely low frequency ballistocardiograph, the recoil 
forces due to cardiac contractions can be calculated by 
referring to the equations in the ballistocardiographic 
literature, substituting the mass of the satellite for that 
of the ballistocardiograph bed and the restoring force 
of the attitude control system for that of the ballisto- 
cardiograph. 
The normal heart pumps approximately 70 ml. of 
blood from each ventricle, 60 to 80 times per minute. 
he pressure varies with each portion of the cardiac 
ycle but the mean pressure in the aorta is about 125 
m Hg. This corresponds to an equivalent rectangular 
eak of blood velocity of about 113 em. per second in 
he aorta, lasting about .2 second per beat. (The 


Mp x Vp = (M) + Mo) Y 
Fie. 2 


actual peak velocity may reach 160 cm./sec. or greater 
for a much shorter period of time.) (6) 

For a 1000 Kg. satellite the displacement due to 
ventricular contraction would be roughly equal to: 


Pe Vo Xm ' 
Ms 
Pia ees 
10° 
where V, = velocity of blood (113 em/sec) 


m, = mass of blood (70 gm.) 

m; = mass of satellite (1000 Kg.) 

¢t = duration of impulse (approx. 0.2 sec) 
This is a trivial displacement except during ultra 
precise work (i.e., photo mapping, etc.) if the operator 
is at the center of gravity of the vehicle, but may 
become more noticeable if the operator is displaced 
from the center of gravity or if the cardiac rate excites a 
system resonance. 

A check on these calculations is available by re- 
computing the empirical normal standards for ampli- 
tude of displacement for low frequency ballistocardio- 
graphs given in the WADC Technical Report 59-593 
“Handbook of Circulation” for a combined body and 
platform weight of 1000 Kg. The actual movement of 
the subject in space, from the tip of the J wave to the 
tip of the M wave averages .72 X 107° cm. This figure 
is quite close to that computed from purely specula- 
tive information. 


Mass Displacement 


It is virtually impossible to specify the center of mass 
or the moment of inertia in three axes of a living 
animal or man. Factors which may affect center of 
mass and movement of inertia and hence influence 
attitude control are, in order of importance: 

a. movements of the arms and feet, 

b. change of the body position within the restraining 

unit, 

c. changes in the distribution of contents within the 

digestive tract, 

d. blood pooling during various physiological proc- 

esses. 

Because of the variation in the position of the axes, the 
center of gravity varies accordingly. For man, the 
position of the center of gravity may, at its best, be 
described as being located in a sphere with a 10 cm 
diameter located in the body at roughly the position of 
the navel. Although there are definite changes in blood 
distribution, these are minor when compared with the 
dynamic influence caused by gross movements of the 
man, shifts of masses caused by digestive processes and 
changes in mass by metabolism (2), (3), (8). 
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Torques or Couples 


The heart and blood vessels form a closed assemblage 
of tubes through which the blood is circulated by the 
action of the heart. The human circulatory system is 
almost, but not quite, bilaterally symmetrical. There is 
a slight dissymmetry of blood flow caused by the bend 
in the aorta which is not compensated for by other 
blood vessels. Since the bend is to the right, the force 
tending to rotate the body would be to the left. 

The heart and circulation can, for the purpose of 
estimating system dynamics, be likened to a simple 
physical circulation model consisting of a pump forcing 
fluid through a network of tubing. Reduced to its 
simplest form, this network is a single elliptical loop 
of tubing with major and minor axes proportionate 
to the dissymmetry of the circulatory system about the 
appropriate axes of the body. 

Let us assume, for the moment, that the momentum 
characteristics of the circulatory system can be simu- 
lated by considering the entire blood flow taking place 
in a loop of tubing with a radius of 5 cm. This distance 
was chosen because cursory measurements of a model 
of the circulatory system indicated that the RMS dis- 
tance of major blood vessels and bends in these blood 
vessels approximated 5 cm. from an arbitrary reference 
point at the tip of the heart. 

The heart of a normal man pumps between 4 and 6 
Kg. of blood per minute; or about 7200 Kg. in a 24 
hour day. (9) At any time, he has the equivalent of a 
fluid flywheel within his body with a mass of approxi- 
mately 5 Kg., a radius of 5 cm., and rotating once a 
minute. 

The angular momentum of a body is defined as the 
produce of its rotary inertia and the angular velocity 
around a given axis. For a flywheel with its weight 
concentrated in the rim, the moment of inertia (J) = 
MR? where M is the mass and F is the radius. Hence 
the angular momentum of the equivalent blood circuit 
ee 


_ © _X_10°)(25)(6.28) 


= MRW 
60 


= 1.31 X 104 gm 
9 = / 
cm? radians/second 


For a 1000 Kg. satellite, of approximately the shape of 
the Mercury vehicle, the moment of inertial along its 
axis 1s that of a right cone of altitude H radius of base R, 
and mass M. This is: 
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For the Mercury vehicle, 1 = 1000 Kg 
R = 100 cm 
H = 300 cnt 


Hence: 
I = 108 X 10* (.3) 


According to the law of conservation of momentu 
the angular momentums of the two rotating naa 
, 
t 


= 3 <x 10? ¢m em? 


must be equal. Hence: 


1.31 X 10‘ gm cm? rad./sec. 3 X 10° gm cm? W © 
Ang. momentum of blood Ang. momentum of va 


W = 4.36 < 10° rad./sec. 


Pa 


or 
2.5 X 10* deg./sec. 


In a 24 hour day, a 1000 Kg. satellite would rotat 
21° 40’ about its axis simply in reaction to the angula 
momentum of the astronaut’s circulatory systen 
Provided the astronaut remained in relatively cor 
stant orientation to the vehicle, as would be the cas 
in most small manned satellites, this torque woul 
affect the attitude of the satellite. Over a long perio 
of time, a significant amount of energy would be re 
quired from the attitude control system to counterac 
this torque. 
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Greek letters and other symbols not available on the 
typewriter should be carefully inserted in ink. Each 
such symbol should be identified unambiguously the 
first time it appears. The distinction between capital 
and lower-case letters should be clearly shown. Avoid 
confusion between zero (0) and the letter O; between 
the numeral (1), the letter 1, and the prime (’); between 
alpha and a, kappa and k, mu and uw, nu and », eta and n. 

The level of subscripts, exponents, subscripts to sub- 
scripts, and exponents in exponents should be clearly 
indicated. 

Complicated exponents and subscripts should be 
avoided when possible to represent by a special symbol. 

Fractions in the body of the text and fractions occur- 
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ring in the numerators or denominators of fractions 
should be written with the solidus. Thus 


cos (x/2b) 
cos (ra/2b) 


is the preferred usage. 

The intended grouping of handwritten formulas can 
be made clear by slight variations in spacing, but this 
procedure is not acceptable in printed formulas. To 
avoid misunderstanding, the order of symbols should 
therefore be carefully considered. Thus 


(a + bx) cos t is preferable to cos t (a + bx) 


In handwritten formulas the size of parentheses, 
brackets and braces can vary more widely than in 
print. Particular attention should therefore be paid to 
the proper use of parentheses, brackets, and braces 
(which should be used in this order). Thus 


{[a + (b + ex)"] cos ky}? 


is required rather than ((a + (b + ex)®") cos ky)?. 
Equations are numbered and referred to in text as 
(15). 


Illustrations 


Drawings should be made with black India ink on 
white paper or tracing cloth, and should be at least 
double the desired size of the cut. Each figure number 
should be marked with soft pencil in the margin or on 
the back of the drawing. The width of the lines of such 
drawings and the size of the lettering must allow for the 
necessary reduction. Reproducible glossy photographs 
are acceptable. However, drawings which are unsuitable 
for reproduction will be returned to the author for re- 
drawing. Legends accompanying the drawings should 
be typewritten on a separate sheet, properly identified. 


Security Clearance 


Authors are responsible for the security clearance by 
an appropriate agency of the material contained in the 
papers. 
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